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Abstract. 

We demonstrate that the processing of the experimental data for the average velocity profiles 
obtained by J. M. Osterlund (www.mesh.kth.se/~jens/zpg/) presented in [ x ] was incorrect. 
Properly processed these data lead to the opposite conclusion: they confirm the Reynolds- 
number-dependent scaling law and disprove the conclusion that the flow in the intermediate 
('overlap') region is Reynolds-number-independent. 
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In a recent issue of the Physics of Fluids Osterlund et al. [ x ] presented results of the pro- 
cessing of experimental data for the average velocity field in the zero-pressure-gradient turbu- 
lent boundary layer, published by the first author of f 1 ] on the Internet 
( Ettp : //www . mesh . kth . x/ ~ j ens/zpg/| ) . The conclusion of [ l ] is very definite (p. 1): 'Con- 
trary to the conclusions of some earlier publications, careful analysis of the data reveals no 
significant Reynolds number dependence for the parameters describing the overlap region 
using the classical logarithmic relation'. 

In the present note we show that the processing of the experimental data in p] is incorrect. 
A correct processing is performed, and the result is the opposite: there is significant Reynolds 
number dependence for the parameters describing the intermediate region and the scaling 
(power) law is valid. We also show where the authors of f 1 ] went wrong. 

The correct processing of the Osterlund data. If a scaling (power) law relates the 
dimensionless velocity U + and dimensionless distance from the wall y + (in the notations of 

t 1 ]): 

(I) U+ = A(y+r (1) 

then in the coordinates lg y + , lg U + the experimental points should lie within experimental 
accuracy along a straight line: lg U + = lg A + a lg y + . Therefore our first step was to 
plot the data from all 70 runs available on the Internet in these coordinates. All 70 runs 
corresponding to different Reg yield the same pattern: in the intermediate region between 
the viscous sublayer and free stream the average velocity distribution consists of two straight 
lines. Three examples are presented in Figure 1 and the Table, all the remaining ones are 
similar and can be found in our detailed report [ 2 ] . Thus the intermediate structure (between 
the viscous sublayer and free stream) consists of two self-similar layers: the scaling law (I) 
is valid for the layer adjacent to the viscous sublayer, and the scaling law 

(II) U+ = B{y+f (2) 

holds in the layer adjacent to the free stream. This means, in particular, that some scaling 
law with Reynolds-number-dependent coefficients is valid in the region (I). The coefficients 
A and a, constant for every run, depend on the Reynolds number which is different for 
different runs, and their variation is substantial (see the Table). 

Once this is established, we investigate whether this scaling law can be represented in 
the form 

U + =( -^lnRe + | y )(y+)Ai (3) 
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obtained by us earlier for flow in pipes. In the case of pipe flows Re was ud/v, where u is 



2 



the mean velocity (bulk flux divided by cross-section area), and d is the diameter. But what 
is Re for the boundary layer? 

The effective Reynolds number Re should have the form Re = UA/u, where U is the free 
stream velocity, v the kinematic viscosity, and A a length scale which we cannot a priori 
identify with the momentum thickness 9, as there is no rationale for such identification. So, 
the basic question is whether one can find for each run a length scale A so that the scaling 
law (3) will be valid for the mean velocity distribution in the first intermediate region. If this 
scale exists then the law (3) is not specific to flows in pipes but may be a general law for 
wall-bounded shear flows at large Reynolds numbers. 

To answer this question we took the values A and a for each run, obtained by standard 
statistical processing of the experimental data in the first intermediate scaling region, and 
then calculated two values In Re\ and In Re 2 by solving two equations suggested by the law 
(3), i 

If the values In Re\ and In Re 2 obtained by solving these two different equations (4) coincide 
within experimental accuracy, then the unique length scale A can be determined so that the 
experimental scaling law in the region (1) coincides with the law (3). Indeed these values are 
close — for all Re e > 10,000, the difference A = (In Re 2 — In Re\)/\n.Re does not exceed 
3%, see the table for the examples in Figure 1 and in [ 2 ] for all runs. This allows one to 
introduce, for large Reynolds numbers, an effective Reynolds number Re, for example by the 
relation 

In Re\ = -(In Re\ + In Re 2 ), or Re = a/ ReiRe 2 , (5) 
2 

the geometric mean of Re\ and Re 2 . This Reynolds number defines the effective length scale 
A, which plays for boundary layer flow the same role as the pipe diameter for flow in pipes. 
Remember that the momentum thickness is calculated by integration of the velocity profile 
obtained experimentally: the calculation of the length scale on the basis of the measured 
velocity profile is not more complicated. Furthermore, the scaling law (3) can be reduced to 
a universal form 

1 , / 2aU + \ , , 
il> = -]nl-= )=lny + (6) 



where a = 2 lr ^ e . This formula gives another way to check the applicability of the Reynolds- 
number-dependent scaling law (3) in the intermediate region (1). Indeed, according to (6), 
in the coordinates In y + ,ip, all experimental points should collapse onto the bisectrix of the 
first quadrant. Figure 2 shows that all data for large Reynolds numbers (Reg > 15, 000, 
24 runs) presented on the Internet collapse onto the bisectrix with accuracy sufficient to 
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give an additional confirmation to the Reynolds- number-dependent scaling law (3). For 
lesser values of Ree a small but systematic parallel shift is observed (see [ 2 ]). One possible 
reason is that in these cases the choice of Re according to (5) may be insufficient because in 
particular at small Reynolds numbers the higher terms of the expansion of the coefficients 
of the scaling law over 1/ln Re could have some influence (see the paper by Radhakrishnan 
Srinivasan [ 3 ]); another possibility is that the relation that was to be proved was used as fact 
in the calculation of the skin friction in equation (6) of f 1 ]. 

Why is the data processing in f 1 ] incorrect? The main argument against the power 
law used by the authors of f 1 ] is the following. They introduce the "diagnostic function" 
(p. 3, right) 

r = £^. (7 ) 

U+ dy+ 1 } 

Their statement, "The function Y should be a constant in a region governed by a power law" 
is correct for a fixed Reynolds number. However, this is not true for the ' diagnostic function 
averaged for KTH data\ which is shown in their Figure 6. 

We invite the reader to look at Figure 1 (the situation with all other runs is the same, 
see our report [ 2 ]). It is clear that for each run T, which is equal to d(lgU + )/dlg(y + ), is 
a constant — look at the straight lines in the first intermediate region! However, (see the 
Table) this constant is different for different runs because the slope of the straight lines is 
Re-dependent! Indeed, the slope in the first region decays with growing Reynolds number. 
It is clear why T obtained by the authors based on averages for KTH data, is decreasing: 
the runs with larger Reynolds number and smaller slopes contribute more at larger y + . 

Due to this incorrect procedure the authors were unable to recognize the Reynolds- 
number-dependent power law in their rather rich database. 

Their determination of the constants of the logarithmic law gives an additional illustration 
of the incorrectness of their procedure. They do it in two ways; it is enough to mention what 
they call the traditional procedure (p. 3, right). They take the data representation in the 
traditional lg y + , U + plane, and calculate the constant k in the logarithmic law 

U + = -ln(y+) + B (8) 

K 

'by fitting a log-law relation for each profile using the following traditional limits of the fit: 
Mi = 50, and M = 0.15'. The result of their fit is shown in Figure 3 (their Figure 5) which 
is far from convincing. 

Apparently embarrassed by this result, the authors extended the upper boundary of the 
viscous sublayer to M 1 = 200 (the traditional value is 30-70), and obtained, again, with huge 
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scatter, the value k = 0.38. For the constant B they obtained the value 4.1. Both values are 
the lowest contenders among the values available in the literature (Nikuradze, k = 0.417, 
5 = 5.89; Monin and Yaglom, k = 0.40, 5 = 5.1; Schlichting, k = 0.40, 5 = 5.5). However, if 
the law (8) is universal, the values of the constants should be identical for all high quality 
experiments! 

Finally, according to the authors of [*], the logarithmic layer extends only over 1/6 of 
the boundary layer thickness. In fact, a better representation is provided by the power law 
in the first region. The upper boundary of this first region is always higher than the upper 
boundary of the logarithmic region presented in f 1 ]. 

In conclusion, careful analysis of Osterlund's data, contrary to the claim in [ x ] , reveals a 
significant Reynolds number dependence for the parameters describing the 'overlap' region 
and confirms the Reynolds-number-dependent scaling law. 
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0.177 


6.24 


12.54 


12.48 


12.51 


0.5 



5 




Figure 1. 
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FIG. 5. The von Karman constant determined by a least-squares fit, with the 
outer limit fixed at 77= 0.15 and the inner limit at M t ; O: KTH, M ■ = 50. • : 
KTH, M f =200. Dashed line: KTH, linear fit, M ( = 50. Solid line: KTH, 
linear fit, M, = 200. □: IIT, M, = 50. ■: IIT, M i — 200. Dotted line: IIT, 
linear fit, Af; = 50. Dash dotted line: IIT, linear fit, M ; = 200. 



Figure 3. 



